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ABSTRACT 

We calculate the emission line spectrum produced by the debris released when a white dwarf (WD) 
is tidally disrupted by an intermediate- mass black hole (IMBH; M ~ 10 2 — 10 5 M@) and we explore the 
possibility of using the emission lines to identify such events and constrain the properties of the IMBH. 
To this end, we adopt and adapt the techniques developed by Strubbe & Quataert to study the optical 
emission lines produced when a main sequence (MS) star is tidally disrupted by a supermassive black 
hole. WDs are tidally disrupted outside of the event horizon of a < 10 5 Mq black hole, which makes 
these tidal disruption events good signposts of IMBHs. We focus on the optical and UV emission lines 
produced when the accretion flare photoionizes the stream of debris that remains unbound during the 
disruption. We find that the spectrum is dominated by lines due to ions of C and O, the strongest 
of which are C iv A1549 at early times and [O in] A5007 at later times. Furthermore, we model the 
profile of the emission lines in the [O in] AA4959, 5007 doublet and find that it is highly asymmetric 
with velocity widths of up to ~ 2500 km s _1 , depending on the properties of the WD- IMBH system 
and the orientation of the observer. Finally, we compare the models with observations of X-ray flares 
and optical emission lines in the cores of globular clusters and propose how future observations can 
test if these features are due to a WD that has been tidally disrupted by an IMBH. 



1. INTRODUCTION 

A star is tidally disrupted when its orbit brings it too 
close to a black hole (BH). If the two bodies are sep- 
arated by less than the tidal disruption radius, Rt ~ 
R+(Mbh /M*) 1 / 3 , where R± is the radius of the star 
and Mbh and M* are the masses of the BH and star, 
respectively, the star is ripped apart because the tidal 
acceleration a c ross i ts diameter exceeds its self gravity. 
iFrank fe Reesl (|1976ft calculated the rate at which a BH 
in the core of a cluster or galaxy would disrupt and in- 
gest main sequence (MS) stars and suggested that the 
associated ultraviolet (UV)/X-ray flare would mark the 
presence of an otherwise undetectable BH. This is pos- 
sible because after the star is disrupted, about half of 
the debris becomes bound to the BH. As this bound 
gas falls back towards the BH, shocks and accretion 
onto the B H produce a luminous flare that peaks in the 



UV/X-ray (|Lacv et al J 1 1981 lR^Jl98ajUlmerlll999l). A 
number of analytic (e.g. , |Phinnevlll989l: Kochanekf 19941 : 
iKhokhlov fe Mehdfl99l fstrubbe fc Quataertil2009ft and 



iBrassart k Luminetl[200l 12010ft . 

Detailed models have also been employed to pre- 
dict the rate at which MS stars are tidally disrupted 
by supermassive BHs. However, uncertainties in the 
supermassive BH mass spectrum, the stellar distribu- 
tion function in galactic bulges, and the relation be- 
tween the two complicate calculations of this rate. For 
example, improved measurements of the mean ratio 
of BH mass to b ulge mass made between the stud- 
ies conducted by iMagorrian k Tremaind (|1999ft and 
IWang k Merrittl (|2004h resulted in a factor of ten in- 
crease in the predicted event rate, raising it to ~ 
10~ 5 yr _1 Mpc~ 3 . This rate is constrained by sys- 
tematic surveys for tidal disruption flares. Histori- 
cally, around ten candidate tidal disruptions were first 
identified by serendipi tous observations of UV/X-ray 
flares in surveys (e.g., iGrupe et al.l [l~995t iBrandt et al.l 
19951: iKomossa fc Baddll99ft iKomossa fc Greinerlll99 " 



Gezari et al.ll2006t iCappelluti et al.ll2009t iMaksvm et al 



nume r ical (e.g.. |Evans fc Kochanek 198 9c ILaguna et aLl 
[19931 IKhokhlov et all 119931: iBogdanovic et all 12004ft 
studies have explored the behavior of the debris after 
a MS star is tidally disrupted and the feasibility of us- 
ing the luminosity and evolution of the flare to search 
for and determine properties of the BH. Others have 
proposed flares and possible detonation of the star due 
to tidal compression as an a dditional observable con- 
sequence of tidal disruptions dCarter fc Luminetl I 1982 L 
1981 [Bicknell fe GingolJll983ULuminet fc Marcldll98i 



Kobavashi et al.l 12004 ). Recent work along these lines 



suggests that the compression will produce shocks capa- 
ble of raising the star's surface te mperature to that of its 
core, resulting in a prompt flare (jGuillochon et al1l2009t 



201CftT Follow up observations of the variability of 
these sources ruled out other explanations and con- 
firmed that the luminosity of the sources declined with 
the same t~ 5 ^ 3 time dependance that models predict 
for the decline in the mas s fallback rate (jHalpern et al.l 
120041 lEsquei et aD 12008ft . Subsequent studies have 
discovered additional tidal disruption candidates and, 
by accounting for the efficiency and sky coverage of 
the surveys, determined that supermassive BHs tidally 
disrupt MS stars at a rate consistent with 10~ 5 — 
10~ 4 yr" 1 galaxy" 1 (iDqnlev et al.1 120021: i Gezari et all 
[20081 [20091: lEsauei et al.l 120081 : ILuo et al.l 12008ft . The 
high energy accretion flares are the brightest and most 
immediate sign of a MS star being tidally disrupted by 
a BH, but they are not the only sign. 
In addition to the UV/X-ray emission, IroosI (UMl) 
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proposed that the unbound portion of the debris can be 
photoionized and produce optical emission lines. The na- 
ture of the optical emissio n lines was explored in d etail 
in the numerical models of iBogdanovic et al.l (l2004f) and 
the an alytic models developed bv IStrubbe fc Quataertl 
(2009, hereafter SQ09). In each of these studies, the 
source of the ionizing photons is the accretion flare. The 
prospect of optical emission from tidal disruptions is ex- 
citing for the following reasons. First, SQ09 showed that 
tidal disruptions are bright enough in the optical band 
that they are potential sources for large transient surveys 
(e.g. The Large Synoptic Survey Telescope, The Palomar 
Transient Factory, and Pan-STARS), which suggests that 
the number of detected candidate events will increase in 
the near future. Second, the emission-line signature of 
such an event helps us identify its nature and distinguish 
it from other, more common transients that can occur 
near the nuclei of galaxies, such as supernovae. 

While flares from tidally disrupted MS stars might be 
a m ore frequent probe of dor mant BHs in galactic nu- 
clei, [Sigurdsson & Rccs (|1997fl pointed out that the BH 
can capture compact objects as well because this region 
is also populated with the products of stellar evolution. 
In particular, white dwarfs (WDs) are also subject to 
the processes that drive MS stars into orbits around the 
central BH with pericenter distances R p < Rt, and can 
therefore be tidally disrupted as well. However, since the 
mass of a WD is comparable to that of a low mass MS 
star and its radius is smaller by a factor of ~ 100, Rt 
is much closer to the BH for a WD than for a MS star. 
Consequently, if M B h > 5 x 10 5 M©, Rt for a WD is in- 
side of the BH's Schwarzschild radius, Rs — 2GMbh/c 2 . 
When Rt < Rs, the WD plunges into the BH horizon 
before being tidally disrupted, making such events un- 
detectable by electromagnetic observations. This makes 
WD tidal disruptions an excellent indicator of the pres- 
ence of an IMBH in the center of a galaxy or globular 
cluster, because if the central BH were more massive than 
~ 10 5 M© we would not observe anything. 

Like the MS case, using WD tidal disruption events to 
probe IMBH requires prior theoretical studies to deter- 
mine how the observed flares are related to the prop- 
erties of the WD-IMBH system. While this scenario 
has not received as much attention as the MS case, 
the early stages of the tidal disruption of a WD in 
an unbound orbit around an IM BH were studied nu- 
merically by iFrolov et al.l (|1994l ). who found that ma- 
terial flowed away from the WD i n supe rsonic jets as 
it was disrupted. iRosswog et al.l (|2009l ) followed the 
evolution of the debris further and found that gas re- 
turns to the BH at a rate that is sufficient to power 
an accretion flare at the Eddington luminosity, LEdd — 
1.3 x 10 41 (Af SH /10 3 M©) ergs" 1 , for up to one year. 
Additionally, since it is supported by degeneracy pres- 
sure, when R p < 0.5Rt the tidal com pression of the WD 
can lead to explosive nuclea r burning (|Luminet fc Pichonl 
Il989t IRosswog et all 120091 ) . The latter authors found 
that, in the most favorable cases, the energy generated 
by thermonuclear burning during the compression phase 
exceeds the binding energy of the WD and is similar to 
that of a Type la supernova. However, if the WD does 
not penetrate too deeply within Rt, the case is simi- 
lar to that of the tidally disrupted MS stars described 
above, in which a portion of the debris is accreted by the 



BH and the rest of the material rema i ns un bound and 
flows away from the BH. iSesana et a l. (2008) calculated 
the light curves of optical and near-UV emission lines 
produced when the unbound debris is photoionized by 
the accretion flare and proposed these lines might be an 
electromagnetic counterpart to the gravitational radia- 
tion emitted as a WD initially on a bound orbit spirals 
into the IMBH. Here, we revise and update these results 
by considering the case of a WD initially on an unbound 
orbit and using a more sophisticated model for the be- 
havior of the debris. 

We apply the analytic description developed by SQ09 
to study optical and UV emission lines from tidally dis- 
rupted MS stars to the case of a WD being tidally dis- 
rupted by an IMBH. We discuss the evolution of the de- 
bris in [§2] drawing primarily from published work. In 
[§3"1 we describe photoionization models and the predicted 
emission lines from these events, including emission-line 
profiles. Finally, we discuss our results in the context of 
resent observations in 54 



2. THE DYNAMICAL MODEL FOR THE DEBRIS 

We adopted the analytic prescriptions of SQ09 to 
model the evolution of the debris after a WD is tidally 
disrupted by an IMBH. Their model describes a MS star 
approaching a supermassive BH on a parabolic orbit with 
pericenter distance R p < Rt- However, since all of the 
relevant quantities depend only on the BH mass, the peri- 
center distance, and the mass and radius of the star, 
the model can also be applied to the case of a WD be- 
ing tidally disrupted by an IMBH. Furthermore, even 
though the WD equation of state is stiffer than that of 
a MS star, the deviations in the tidal disruption process 
and the behavior of the debris that arise from this differ- 
ence are small, so the approximate formulae developed in 
SQ09 give a reasonable description of the WD case. The 
model assumes that half of the material becomes bound 
to the BH and forms an accretion disk after the WD is 
disrupted, and that the other half of the debris is ejected 
on hyperbolic orbits. The details of these components in 
the SQ09 model and the adaptations we have made are 
described below. For a full derivation of the model, we 
refer the reader to SQ09. 

The evolution of the debris was modeled for a range of 
initial conditions. To explore how the optical emission 
lines are affected by the mass of the IMBH, we com- 
puted models with M B h = 10 2 , 10 3 , and 10 4 M©. Even 
though a WD of the mass used in the models would be 
tidally disrupted outside of the event horizon of a 10 5 M© 
IMBH, we do not consider this case because the tidal dis- 
ruption radius is within the last stable orbit around the 
IMBH. We also constructed models for two different val- 
ues of the pericenter distance, one set with R p = Rt 
and, for co mparison, a sec o nd se t with R p = 0.3i?r, 
even though IRosswog etail (|2009j _ report that the WD 
would explode in the latter case. ITablc II lists the initial 
conditions for each model we considered. In each of the 
models, we set M\yd = 0.55 M© because the mass distri- 
bution of WDs observed in the S loan Digital Sky S urvey 
is strongly peaked at this value (|Madej et al.ll200ll ). We 
calculated the WD's radius to b e Rwd = 8 - 6 x 10 8 cm 
using the mass-radius relation of iNauenbergj (|1972[) . 

2.1. The Accretion Disk 
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As the bound material falls back to pericenter, it 
shocks on itself and the orbits of the debris circularize to 
form an accretion disk. In the model, the disk extends 
from the last stable orbit, RlsOi out to approximately 
2R p . Like SQ09, we considered both Schwarzschild BHs 
with Rlso — 3Rs and maximally spinning BHs with 
Rlso — Rs/%- The rate at which the bound material 
falls back towards the BH is set by the time it takes 
the debris to return to pericenter and the mass of the 
WD, and was computed using equations (1) and (2) 
of SQ09. The mass fallback rate decreases as t~ 5 / 3 , 
and at early times it can exceed the Eddington mass 
accretion rate, MEdd = WLEdd/c 2 . We capped the 
rate at which material is accreted by the BH, M, at 
MEdd and, following SQ09, we assumed that the excess 
mass is blown off in a radiation driven wind. Because 
the nature of these outflows is uncertain, we did not 
consider emission from the outflows or how the outflows 
might affect the emission lines produced in the unbound 
material. SQ09 suggested that the material in the out- 
flows can also be photoionized and produce additional 
line luminosity, making the predi ctions of this model 
lower limits. In a followup paper, iStrubbe &: Quataertl 
(2010) consider the photoionization of the out-flowing 
gas that is outside the photosphere of the outflow 
and calculate the resulting absorption-line spectrum. 
Furthermore, the outflows could alter the radiation 
field that photoionizes the unbound material and, thus, 
change the ionization state of the unbound debris. Since 
these effects were not accounted for in our model, the 
emission line luminosities predicted during the super- 
Eddington phase are uncertain. The outflows subside 
after a time tedd, given in equation (3) of SQ09 as tEdd ~ 
0.1 (A W10 6 M w ) 2 / 5 (R P /3R S ) 6 / 5 (M*/ M Q ) 3 / 5 (R*/R Q ) 
yr. In ITablc 11 we list the value of tEdd for each of our 
models. 



2.2. The Unbound Debris 

As the star approaches the BH, it is tidally stretched 
along the direction of its motion and then spun up, re- 
sulting in a spread in the specific energy of the debris 
(|Reeslll988l ). SQ09 calculated how this range of specific 
energies leads to an expanding, curved wedge of unbound 
debris spanning a range of azimuthal angles and radial 
distances from the BH (see |Figurc 1\ . The accretion flare 
illuminates the inside face of this wedge, and the gas 
along this edge reprocesses the radiation and produces 
the optical emission lines we are interested in. SQ09 
modeled the extent of this face in terms of i? m ax, A</>, 
and Ai, the maximum radial distance from the BH, the 
azimuthal dispersion of the debris in the orbital plane, 
and the angular dispersion of the debris perpendicular 
to the orbital plane, respectively. We adopt their equa- 
tions for the angular dispersion of the unbound material, 
however, we make a minor modification to calculate the 
distance between the BH and material along the debris 
arc. Since we considered several points along the arc, not 
just the most distant material, we calculated the distance 
to the BH at any point along the arc with 

R(4>) -Uptcot^V (1) 



where v p is the pericenter velocity of the WD, t is the 
time since pericenter passage, and <p is the angle in the 
orbital plane between pericenter and the point in the 
debris cloud being considered and is in the range (jr — 
A0) < <fi < ir. The WD is initially on a parabolic orbit, 

so v v = {2GM BH /Rp) 1/2 - SQ09 calculated R max in a 
similar manner and our equation 01 is c onsistent with 
the derivation of lEvans I fc Kocha nck (1989j) for the values 
of <fi considered here. In ITableTl we list A(f> and w m ax/c, 
the maximum velocity in the unbound debris cloud as 
a fraction of the speed of light, for each set of initial 
conditions. In addition to R and <fi, we define the £-curve, 
which begins at the IMBH and lies along the illuminated 
face of the debris cloud. 

Furthermore, because we computed photoionization 
models at several points along the inside face of the cloud 
of unbound debris, we also had to deviate from SQ09's 
model whe n we calcu l ated t he density of the unbound 
debris, n. iKochanekl (|1994h found that the density in 
the stream of unbound material was not uniform and de- 
creased with distance from the BH asnoc R~ 3 because 
the material is undergoing free expansion. We used this 
density profile in our model and normalized it using the 
fact that the mass of the wedge of unbound material was 
half the mass of the WD. 

3. PREDICTED EMISSION 

Predicting the emission lines we will observe when a 
WD is tidally disrupted by an IMBH requires a model for 
the emission from both the bound and unbound portions 
of the debris. To calculate the luminosity and spectrum 
of the accretion flare generated by the disk of bound gas, 
SQ09 assumed a "slim disk" model (|Abramowicz et aLl 
1 19881 ) and derived an expression for the effective tem- 
peratures as a function of radius in the disk (see their 
equation (19)). Then, using the range of temperatures 
throughout the disk, they modeled the spectrum as a 
multicolor blackbody. We considered this spectral energy 
distribution (SED) for the disk as well as two others. We 
also used a SED that consisted of the multicolor black- 
body and an additional X-ray component from 0.1 — 100 
keV with a power law of the form Je oc E~ 1 , where /e 
is the flux density per unit energy and E is the photon 
energy. We normalized the intensity of the X-ray power 
law using the intensity of the multicolor blackbody at 
25 00 A and the ctg x parameteiQ] based on the findings 
of iSteffen et~aTl (|2006t ). The third SED we considered 
was an empirically constructed, multi-component model 
for the SED of an active galactic nucleus ( AGN) that was 
first used inlKorista et al.l (119971) and explored in greater 
detail bv lCasebeer et al.1 {2006). This model is param- 
eterized by the temperature of the peak of the model's 
thermal component, aox, and the slope of the the X-ray 
component. We set these parameters to the temperature 
of the peak of the multicolor blackbody, —1.4, and —1, 
respectively. The accretion disk-IMBH system is similar 
to an AGN, so we used these three SEDs to sample the 
range of UV/X-ray and UV/infrared ratios observed in 

1 The optical-to-X-ray spectral index is the exponent of a 
power-law connecting the monochromatic luminosity densities at 
2500 A and 2 keV. If L„ oc u~ a ox _ then a x = 

- [log L v (2500 A) - logZ„(2 keV)] / [logv(2500 A) - Iogf(2 keV)] = 
1 + 0.384 [log (vL„) 25oo k - ]og(vL u ) 2 kcV ] • 
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AGN. This is necessary because heating of the free elec- 
trons in the unbound material by infrared radiation can 
change its ionization structure and irradiating this ma- 
terial with X-rays can lead to significant heating per ion- 
ization, which would affect the luminosities of collision- 
ally excited lines. In all three cases, the time evolution 
of the SED is governed by the declining mass accretion 
rate, which decreases the bolometric luminosity of the 
SED and causes its peak to shift to longer wavelengths 
as time goes on. |Figure 2] shows each SED and its evolu- 
tion. These SEDs were used as one of the inputs to the 
photoionization models described in the next section. 

3.1. Photoionization Models 

The unbound material is assumed to radiate because 
it is illuminated by the accretion flare described above. 
Predicting the emission lines generated by this material 
requires photoionization models to determine how the 
gas reprocesses the incident radiation. We performed 
the photoionization calcul ations with version 08.00 of 
Cloudy, last described by iFerland et al.l (|1998t ). Since 
the debris cloud consists of material from a tidally dis- 
rupted WD, its composition is different from that of a 
MS star. The composition of WDs is predicted by the 
yields of nuclear reactio ns and confi rmed with asteroseis - 
mology (|Metcalfdl2003l ). Following iMadei et al.l (|2004D . 
we assumed mass fractions of 67% O, 32% C, and 1% 
He and all other elements with their relative abundances 
scaled from the solar values. Hydrogen is assumed to 
make up 0.001% of the WD's mass. For eac h set of ini- 
tial conditions (i.e., models A-F in ITableTj) . we calcu- 
lated the state of the accretion flare SED and the wedge 
of unbound debris at a series of times after the tidal dis- 
ruption. For each time step, we split the wedge into 
six azimuthal segments. With six segments, we were 
able smoothly interpolate quantities along the I direc- 
tion. Next, we calculated the density of the gas and the 
intensity of the ionizing radiation at the center of each 
segment. Several considerations went into the calcula- 
tion of the latter. First, each segment is at a different 
distance from the accretion disk and has a correspond- 
ingly different geometric attenuation factor. Next, since 
the segments lie along an arc, the ionizing radiation's an- 
gle of incidence is also different for each segment. Finally, 
in addition to changes in the intensity of the incident ra- 
diation, the shape of the SED incident on each segment 
also varies. The light travel time, ti t , to the unbound de- 
bris furthest from the BH is comparable to the timescale 
on which the accretion flare evolves, so the SED of the 
radiation incident on these portions of the wedge was 
emitted at t — tu, when the peak was at a higher en- 
ergy. Each of these considerations has a significant effect 
on our models. Geometric attenuation and light travel 
time dominate for more distant regions, but the material 
nearest the disk is illuminated at nearly glancing inci- 
dence, significantly reducing the intensity of the incident 
radiation despite the proximity to the disk. With these 
parameters, we are able to calculate how the unbound 
debris responds to the ionizing radiation emitted by the 
accretion disk. 

The ionization state of the gas is convention- 
ally described by the ionization parameter Uh = 
Qh/(^it-hR 2 c), where Qh is the emission rate of pho- 
tons with energy greater than 1 Ry, nn is the hydro- 



gen number density, and R is the distance to the ion- 
izing source. However, since hydrogen is depleted in 
the material considered here, it is advantageous to use 
a more abundant reference element when calculating 
the ionization parameter. We chose to use oxygen be- 
cause its first ionization potential is very near that of 
hydrogen so the ratio of the two ionization parameters 
is Uo/Uh = Xh/Xo, where Xo and Xh are the oxy- 
gen and hydrogen abundances by number, respectively. 
Thus, Uo is the ratio of the density of photons capable 
of singl y ionizing oxygen to the total density of oxygen 
atoms. |Figure 3 shows how Uo, the total number den- 
sity of all elements, and the radial distance to the IMBH 
change along the illuminated face of the debris cloud, £, 
for model C at 600 and 1200 days after tidal disruption. 
Since the density drops faster with distance from the ac- 
cretion disk than the flux of ionizing photons does, Uo 
increases along the ^-curve. While the fallback rate is 
super-Eddington, the value of Uo increases linearly with 
time at all points along the £-curve because the luminos- 
ity of the ionizing source remains constant at LiEdd, while 
the density of the expanding debris cloud decreases. We 
find numerically that Uo declines approximately as t~ 0/2 
after the super-Eddington accretion phase endtQ. 
The fractions of the first four ions of oxygen and carbon 

along h from a represen tative photoionization calculation 
are shown in |Figure 4) In nearly all of our photoioniza- 
tion models with R p = Rt (i.e., models A,C,E, and F), 
the dominant ion of oxygen is O n and the dominant ion 
of carbon is C n. The exceptions are the azimuthal seg- 
ment nearest the IMBH at times less than 300 days after 
tidal disruption, which consists mostly of neutral oxygen 
and carbon, and the segment farthest from the IMBH 
at times greater than 2500 days after tidal disruption, 
whose dominant ion are O in and C m. As the den- 
sity and ionization parameter change, the relative thick- 
ness of the zone dominated by each ion of oxygen does 
not change significantly between 300 and 2500 days after 
tidal disruption, so the relative abundances of the ions of 
oxygen do not vary significantly in these models. In the 
case of carbon, however, the abundance of C in relative 
to C iv decreases during the first 600 days after tidal 
disruption, before leveling off. During this time some of 
the C in is photoionized into C iv. For the models with 
R p = Q.3Rt (i.e., models B and D), O n and C n are 
the dominant ions of oxygen and carbon, respectively, 
throughout the debris tail for the first 300 days. In these 
models, the density of the unbound debris is much lower, 
so Uo is much higher. Therefore, at later times, higher 
ionization species (O iv, O v, C iv, and C v) dominate 
in different portions of the cloud. 

In the photoionization models, we make the simplify- 
ing assumption that the density of the debris is constant 
along h, the path the ionizi ng radiation follows into the 
debris cloud (see |Figure lj ). This assumption is valid 
because the thickness of the ionized layer is very small 



2 This agrees with a simple analytic estimate, Uo is the ratio 
of the density of ionizing photons to the density of oxygen atoms 
so Uo oc L c i/(T c iR 2 n) oc 4 -1 / 4 where and are the bolomet- 
ric luminosity and maximum temperature in the accretion disk, 
respectively. The time dependance is slightly steeper than in our 
models because this estimate does not adequately account for the 
broad peak in the SEDs. 



EM Emission From Tidally Disrupted WD 



5 



compared to the distance to the BH, so the density in this 
skin only deviates from the surface value by 0.1% in the 
most extreme cases. In addition, we assumed that the 
ionized gas remains in photoionization equilibrium dur- 
ing the 4000 days probed by our models. This requires 
that the density of the gas and flux of ionizing radia- 
tion change on timescales greater than the recombination 
time t rec ~ (n e a rec ) _1 of the ions of interest. Since the 
cloud consists mostly of carbon and oxygen, there are 
many electrons per ion, which keeps the electron density 
reasonably large even at late times when the density of 
the debris cloud is low. Furthermore, the range of condi- 
tions encountered in the models yielded O iv, O III, O u, 
C iv, and C in recombination coefficients that, in combi- 
nation with the electron density, kept the recombination 
time well below the timescales on which the conditions 
in the wedge of unbound debris changed (the ranges of 
recombination rates are 6.7 x 10 -7 — 1.9 s^ 1 , 1.3xl0~ 7 — 
0.52 s _1 , 4.4xl0~ 8 -0.15 s _1 , 1.3 x 10~ 7 -0.97 s" 1 , and 
3.8 x 10~ 7 — 0.98 s _1 for the above ions at the lowest and 
highest electron densities, respectively). Therefore, the 
cloud can remain in photoionization equilibrium for over 
ten years. 

3.2. Line Profiles 

Combining the dynamical and photoionization models 
allows us to calculate the observed profile of the emission 
lines in the [O in] AA4959, 5007 doublet. The dynami- 
cal model we have adopted gives us an analytic expres- 
sion for the velocity of the material along the illumi- 
nated edge of the debris cloud (£) in the orbital plane, 
v r{<P) ~ Vp cot(</>/2). This is the dominant component 
of the debris cloud's velocity, so we neglected the motion 
of the material perpendicular to the orbital plane, which 
is smaller by a factor of Rwd/Rp (SQ09). The veloc- 
ity was projected onto the observer's line of sight, which 
is defined by io and 9q, the angles between the line of 
sight and the orbital plane and the pericenter direction, 
respectively. |Figure l" shows the io = $o = axis. We 
interpolated the [O mj A5007 flux computed in our pho- 
toionization models along the £-curve to determine the 
restframe flux at each velocity along I. We took the flux 
at 4959 A to be one third of this value. We considered 
Doppler boosting of this emissivity, but the effect is in- 
significant at most points along the arc and produces a 
maximal change of 1.2%. Finally, at each point, we took 
into account the relativistic Doppler shift, the gravita- 
tional redshift, and local broadening. We found that the 
gravitational redshift was negligible in all of our models. 
The local broadening was assumed to have a Gaussian 
profile and included terms for the thermal motions of 
the gas, the range of velocities within each bin along the 
£-curve, and the range of velocities in the ionized skin, 
along h. Since the conditions along the ^-curve varied 
significantly, each source of local broadening considered 
in our models dominates in a different portion of the 
illuminated face of the cloud at any given time. Fur- 
thermore, the conditions in the unbound debris change 
substantially over time, so the dominant source of local 
broadening in a given portion of the debris cloud also 
changes over time. In the next section, we discuss the 
results of the photoionization and emission-line profile 
calculations. 



3.3. Results 

Our photoionization calculations predict an emission- 
line spectrum that is dominated by lines from carbon and 
oxygen. The light curves of some of the strongest emis- 
sion lines are shown in |Figure 5] and the luminosities of 
th e six stro ngest features at two different tim es are lis ted 
in ITable 21 The values given in Figure 5| and ITable 21 are 
from calculations preformed with the multicolor black- 
body plus X-ray power law SED, but calculations with 
the multicolor blackbody SED produced similar results 
because the ionization state of the gas is largely deter- 
mined by the flux of UV photons, which is nearly iden- 
tical in these two SEDs. In the case of the multicompo- 
nent, empirical AGN SED, the emission-line luminosities 
we calculated were consistent with those calculated using 
the other two SEDs, only varying by a factor of ~ 2 in the 
most extreme cases. The AGN SED has much stronger 
optical and near-IR continuum emission, which allows 
the continuum to outshine the emission lines, greatly re- 
ducing their equivalent widths. However, these results 
indicate that heating of the free electrons in the unbound 
material by IR radiation from the IMBH-accretion disk 
system does not have a strong effect on the emission line 
luminosities. Thus the main difference between the re- 
sults from the AGN SED and the other two cases is that 
the optical continuum is higher and equivalent widths of 
the optical emission lines are lower, even though their 
luminosities are approximately the same. Our discussion 
will focus on the spectra predicted by calculations using 
the other two SEDs, which are nearly identical in the 
optical and UV bands. 

For each model with R p = Rt (i.e., models A, C, E, 
and F), the UV lines C iv A 1549 and C in A977 have the 
highest luminosity early on. The luminosities of these 
and other permitted lines decline over time. The lumi- 
nosity of the C in A977 line declines faster than the C iv 
A1549 line because of the reduction in the amount of 
C in relative to C iv discussed in §3.1. After 100 days, 
the [O in] A5007, [O in] A4363, and [O n] A7325 lines 
have the largest equivalent widths. These are forbidden 
transitions, which are collisionally excited and then de- 
excite radiatively as long as the electron density is below 
the threshold for collisional de-excitation (i.e., the crit- 
ical density for the transition n cr , listed in Table [5] for 
reference). In the early phases of the evolution of the 
debris (the super-Eddington phase), the luminosity of 
these forbidden lines rises with time and the lines reach 
their maximum luminosities at 600, 250, and 300 days, 
respectively, when the outer tip of the ionized debris tail 
is at the corresponding critical density. After the peak, 
the line luminosities begin to decline, like the permitted 
lines. To understand this trend, we note that, at den- 
sities below the critical density of a transition the line 
emissivity per unit volume increases linearly with den- 
sity and peaks at the critical density. As the density in- 
creases above the critical density, the line emissivity per 
unit volume drops by a factor of several (with the deple- 
tion of the upper level population by collisions) and then 
levels off (as the upper level is re-populated by collisions 
and radiative decays from higher levels) Thus, the ini- 

3 We verified this behavior, by solving the population equations 
for the [O III] A5007 transition and then combining the result- 



6 



Clausen & Eracleous 



tial rise in the luminosity of the forbidden lines is the 
result of the expansion of the debris while the density is 
above the critical density of the transition. This expan- 
sion occurs during the super-Eddington phase, when the 
ionizing luminosity is constant. During this expansion, 
the peak luminosity in the light curves corresponds to 
the time when the density at the outer tip of the debris 
tail reaches the critical density. After that time the den- 
sity of an increasing portion of the debris drops below 
the critical density (causing a reduction in the emissivity 
per unit volume) and the total luminosity declines. We 
can quantify the rate of the initial rise in a forbidden 
emission line's luminosity if it occurs while the ionizing 
luminosity is constant and the density in the debris cloud 
is above the critical density. As the debris expands, the 
column density of the ionized skin remains constant be- 
cause the decline in the flux of ionizing photons resulting 
from the increasing distance from the IMBH is balanced 
by the increase in the cloud's surface area. This column 
sustains a constant emission line flux at the surface of 
the cloud, so the rise in the forbidden emission line lumi- 
nosity is driven by the cloud's expanding emitting area. 
This area grows as i 2 , and this evolution can be seen in 
Figures [5] and HI 

We show the [O m] A50 07 and C iv A 1549 light curves 
for each of the models in | Figure 6] These curves come 
from calculations preformed with the multicolor black- 
body plus X-ray power law SED and we note that cal- 
culations that exclude the X-ray component of the inci- 
dent radiation predict luminosities within a few percent 
of those shown. We calculated the luminosities by mul- 
tiplying the emergent flux computed by Cloudy in each 
azimuthal segment by that segment's emitting area and 
then summing the luminosities across the entire illumi- 
nated face. The light curves indicate how the evolution 
of the [O in] A5007 line depends on the properties of the 
WD-IMBH system and how its evolution can be used to 
identify such a tidal disruption event. First, the lumi- 
nosity of the line at any given time increases as the mass 
of the IMBH decreases. This is because both the initial, 
maximal mass accretion rate and the emitting area of 
the cloud of unbound debris decrease as the mass of the 
IMBH increases. Next, the line reaches its peak lumi- 
nosity at earlier times for smaller values of R p /Rt (i-e., 
models B and D). This is primarily because the spread in 
the debris is larger for smaller R p and the critical density 
of the transition is reached earlier. The steeper decline 
of the [O in] A5007 luminosity in models B and D is also 
a result of the lower density in these debris tails. The 
less dense debris becomes highly ionized, reducing the 
amount of O in in the cloud. In any case, the decline 
of the line luminosity after it reaches its maximum value 
is less step than the t~ 5 ^ 3 decline of the mass fallback 
rate. Finally, if the IMBH is spinning, the inner radius 
of the accretion disk is much closer to the IMBH and 
therefore has a maximum temperature that is 2.7 times 
higher than that of a disk around a Schwarzschild BH. 
This leads to an increase in the flux of ionizing photons 
and a corresponding overall increase in the luminosity of 

ing level populations with the luminosity per unit density of ions 
in the lower level of the transition. To carry out this exercise, 
we assumed a 5-level i on an d followed the methodology in §3 of 
IQsterbrock fc Ferland] (2006ft . 



the [O in] A5007 line. 

The emission-line profile of the [O m] AA4959, 5007 
doublet depends on the properties of the WD-IMBH 
system. Qualitatively, the full width at half maximum 
(FWHM) of each line is relatively small shortly after the 
tidal disruption. At these times, the debris cloud is above 
the critical density and the line profile is dominated by 
emission from the outer tip of the debris cloud where 
the emitting area is maximal. The spread of velocities 
in this small segment of the cloud is comparatively nar- 
row, resulting in an emission line with a small FWHM. 
This material has the greatest velocity with respect to 
the IMBH, so the peak of the emission line will be off- 
set from its rest value. Eventually, the density of the 
material at the outer tip of the debris tail drops to n cr 
and the emission line reaches its maximum luminosity. 
As time goes on, material at smaller £, that is moving 
at slower velocities relative to the IMBH, reaches n cr for 
the [O in] AA4959, 5007 doublet and the emissivity per 
unit volume increases. At the same time, the increas- 
ing volume of material in the more distant regions of the 
cloud partially compensates for the decreasing emissivity 
per unit volume of this material. As a result, material 
at a large range of velocities produces the [O in] AA4959, 
5007 doublet at comparable luminosities and the FWHM 
of each line increases. The broadening is asymmetric be- 
cause the velocity of the gas along the illuminated edge 
increases monotonically with distance from the IMBH. 
The orientation of the observer determines the degree 
of asymmetry and whether the wing blueward or red- 
ward of the peak is broader. Furthermore, the position 
of the peak of each line also shifts as the debris cloud 
expands and the region of maximum emissivity moves 
toward the IMBH along the £-curve. The emissivity is 
largest when and as the cloud expands, regions 

near the IMBH traveling at lower velocities relative to 
the IMBH reach this density and become brighter at 5007 
A. Later, when most of the unbound debris is below the 
critical density, the brightest region can shift away from 
the IMBH towards a region with a larger volume, and 
in some cases a larger fraction of O in. Again, whether 
the peak shifts towards the red or the blue is determined 
by the orientation of the observer. The range of veloci- 
ties along the inner edge of the debris cloud is broader for 
higher mass BHsQ , but since the FWHM of observed line 
profiles depends on the velocity along the line of sight, 
broader lines do not necessarily mean more massive BHs. 

To illustrate these points, Figures [7] and [8] show the 
time evolution of the [O in] A5007 line for model C, as 
observed from 8 = —5° and 6 = 30°, respectively, with 
i{) = in both figures. In the former, the observer is 
positioned such that the spread in line of sight velocities 
along the brightest part of the cloud is quite large. In the 
latter, the velocity of the material is mostly transverse 
to the line of sight, so the range of observed velocities in 
the brightest part of the cloud is relatively small. These 
figures do not include the [O in] A4959 line, so the dis- 
tinct "components" in the line profile result from the 

4 SQ09 showed that the maximum velocity of the unbound ma- 

1/2 1/3 

terial is u max oc v p R p . For tidal disruption, R p < Rt oc M bh . 

Since v p oc (Mbh /Rp) 1 ^ 2 > the maximum velocity of the unbound 

material increases as Af!/J 
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distribution of [O m] A5007 emitted flux as a function of 
projected velocity. Unfortunately, since changes in the 
position of the observer result in such drastic changes in 
the line profile, it is difficult to extract detailed infor- 
mation about the nature of the WD-IMBH system from 
it. However, as we discuss in the following section, the 
general properties of the emission-line profiles are useful 
in identifying candidate WD tidal disruptions. 

4. DISCUSSION 

We have presented calculations of the line emission pro- 
duced when the unbound debris of a WD that has been 
tidally disrupted by an IMBH is illuminated by the asso- 
ciated accretion flare. We found that the two strongest 
emission lines, by far, are C iv A1549 at early times and 
[O in] A5007 at late times. Furthermore, the emission 
lines formed in the illuminated face of the debris cloud 
will have broad, asymmetric profiles with velocity widths 
of several hundred to a few thousand km s _1 . Our results 
suggest that UV and optical spectra whose strongest fea- 
tures are broad C iv A1549 and [O m] A5007 emission 
lines can serve to identify WD tidal disruption events, 
when they are detected in addition to an X-ray flare in 
the cente r of a globular cluster or dwarf elliptical galaxy. 
|Figurc 6] shows how the characteristic timescales and lu- 
minosities of these emission lines scale with R p and the 
mass of the IMBH. The luminosities of the C iv A1549 
and [O in] A5007 lines decline more slowly if R p = Rt 
than they do if R p = O.SRt- For a constant R p , the 
shapes of the emission-line light curves do not change 
much when the mass of the IMBH is changed, however, 
systems with lower mass IMBHs produce brighter C iv 
A1549 and [O in] A5007 emission. The emission line lu- 
minosities are sensitive to the density of the photoion- 
ized debris, and since this density is only weakly de- 
pendent on Mbh, the line luminosities evolve on simi- 
lar timescales despite a factor of 100 increase in Mbh- 
In these models n = no/R 3 , where the normalization 

factor n = M%% M B % (8/x)- 1 {R P /3R T ) 1/2 {tan[(7r - 
A0)/2][sinAi (log(sm(A<£/2)) + 6)]" 1 } ensures that the 
mass of the unbound debris tail is 0.5Mwd- Here we 
have used the mean molecular mass \x. Numerical fits to 
the quantity in the curly brackets show that it increases 
with the mass of the IMBH as M^ 4 when R P /R T is held 
constant. Altogether, we find that the density of the un- 
bound debris is only weakly dependent on the mass of 
the IMBH, n cx Mjjjy , so it is difficult to constrain Mbh 
with the [O in] A5007 light curve. 

Our work drew heavily on the SQ09 model for the 
tidal disruption of a MS star by an supermassive BH, so 
we summarize the uncertainties in the model described 
there. First, since the accretion disk becomes geomet- 
rically thin after M drops below the Eddington rate, it 
becomes susceptible to viscous instabilities that impede 
the steady flow of material through the disk and this 
could result in a deviation from the smooth i -5 / 3 depen- 
dance of the ionizing radiation. Material will accumulate 
in the disk rather than accreting onto the BH when disk's 
viscous timescale becomes comparable to the time since 
tidal disruption. This occurs at t ~ 8 — 12 yr for the disks 
in our models. Furthermore, the model is no longer valid 
when the mass accretion rate drops below ~ O.OlMEdd 
and the disk becomes radiative ly inefficient. The final 



time steps in our models B-F are nearing this limit, and 
this change will lead to a reduction in the number of 
UV photons not accounted for in this model. Since the 
photoionization cross-sections for O and C decline signif- 
icantly from UV to X-ray energies, the ionization state 
of the unbound debris depends sensitively on the flux of 
UV photons whose energies are near the ionization po- 
tentials of C and O. Next, a significant amount of initially 
bound material can be driven away from the IMBH in the 
super-Eddington outflows. The impact that the outflows 
have on emission from the unbound material is uncertain. 
SQ09 pointed out that the outflows could also be pho- 
toionized by the accretion flare and generate additional 
line luminosity. Additionally, since we did not account 
for how the outflows effect the radiation that illuminates 
the unbound debris, the emission line luminosities pre- 
dicted for t < ^Edd are uncertain. 

The modifications we made to the SQ09 model allowed 
us to test how some of their simplifications affect the 
predictions of the model. In their photoionization calcu- 
lations, they simplify the geometry of the wedge of un- 
bound debris by approximating it as a cloud separated 
from the ionizing source by i? max (i-e. the maximum of 
equation ([T]) or R(tt — A<fi)). Furthermore, they assume 
that the density of the unbound material is uniform and 
evaluate it by placing half the mass of the star in a re- 
gion of volume i?^ AR Ai A<fi. In our calculations, on 
the other hand, we have accounted for the changing ion- 
ization state of the material along the illuminated edge 
of the cloud by allowing the density and intensity of the 
ionizing radiation to vary along this edge, as discussed in 
§3.1. Both methods produce similar emission-line light 
curves, with the SQ09 model predicting slightly higher 
luminosities for most lines. This was true for both the 
WD-IMBH case and the MS-supermassive BH case. Our 
models also show that the choice of SED for the accretion 
flare does not have a strong affect on the emission line 
luminosities. We conclude, therefore, that their simpler 
method is sufficient to predict the luminosities and light 
curves of optical emission lines from tidally disrupted 
stars. 

One additional source of uncertainty in our models is 
th e composition of the unbo und material. The models 
of iLuminet fe Pichonl (|1989f) and iRosswog etaH (|2009f ) 
found that compression of the WD as it passes through 
pericenter resulted in nuclear burning. The amount of 
energy released by these reactions and how significantly 
they alter the composition of the debris depends on the 
initial composition of the WD, Mbh, and R p . To ex- 
plore how a change in the abundances brought about by 
nuclear burning would influence the emission lines emit- 
ted by the unbound debris, we also comp uted photoion- 
ization models using abundances given in Rosswog et al. 
(2009) for a 1.2 M Q WD (R WD = 2.7 x 10 s cm) that is 
tidally disrupted by a 500 M© IMBH with R p = R T /3.2. 
For the first year, the brightest lines are the UV lines C iv 
A1549, Si in' A1888, and O vi A1035. Later however, the 
emission-line spectrum is dominated by forbidden iron 
lines rather than carbon and oxygen lines. The most lu- 
minous lines include [Fe x] A6373, [Fe vn] A6807, [Si vm] 
A1446, [Fe vn] A3759, [Fe v] A3892, and [Fe vn JA5271. 
While, this scenario represents an extreme case in which 
the energy released by nuclear burning is greater than 
the binding energy of the WD, the presence of these iron 
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and silicon lines in a spectrum that is otherwise simi- 
lar to that described above could be evidence of nuclear 
burning during the tidal disruption. 

The next question to consider is how often we should 
expect to see emission lines fr o m tid ally disrupted 
WDs. iRamirez-Ruiz fc Rosswoj <|2009f ) estimate that 
the rate of LEdd flares from MS stars being tidally dis- 
rupted by 10 3 - 10 4 M Q IMBHs in globular clusters is 
~ 4000 yr -1 Gpc~ 3 . Following these authors, we as- 
sume all globular clusters have an IMBH in their cen- 
ter and a dopt a globular cluster space density of uqc ~ 
4 Mpc" 3 ([Brodie fe Straderll2006l ). We take the WD tidal 
disruption rate of 10~ 8 yr -1 (globular cluster) -1 from 
iSigurdsson fc Ree"sl ([19971 ) and arrive at a WD tidal dis- 
ruption rate of ~ 40 yr _1 Gpc -3 , 100 times lower than 
the MS rate. This rate is a very optimistic estimate be- 
cause it is unlikely that every globular cluster has an 
IMBH in its center and the actual rate could be much 
lower. The likelihood of observing such an event is fur- 
ther reduced by its relatively short lifetime of ~ 10 yr. 

Des pite the meager rate and limited observing win- 
dow, llrwin et al.l (|2010f ) proposed the tidal disruption of 
a WD by an IMBH as an explanation for two recent ob- 
servations of X-ray a nd optical emission in the centers 
of globular clusters, llrwin et"a"l . (2010) report on the 
ultraluminous X-ray source CXOJ033831.8-352604. The 
source has a 0.3 — 10 keV luminosity of ~ 2 x 10 39 erg s _1 
as well as [O m] A5007 and [N n] A6583 emission lines 
with luminosities of few x 10 36 erg s _1 and FWHMs of 
140 km s _1 . The authors discuss a second, similar source 
in the center of a globular clu ster associated with NGC 
4472. Maccaro ne et al.l (|2007| ) conclude that the source's 
X-ray luminosity of 4 x 10 39 erg s _1 and its variability 
can only be explained b y the presence of a BH. In follow 
up optical observations. IZepf et al.l ([20081 ) found that the 
[O in] A5007 line luminosity was 1.4 x 10 37 erg s _1 and 
that the line's FWHM was 1500 km s -1 . 

In both cases, the X-ray luminosity of the source is 
consistent with the models that assumed a multicolor 
blackbody plus X-ray power law SED for the accretion 
flare. However, in all of our models, the luminosity of the 
[N n] A6583 line is two orders of magnitude lower than 
that of the [O in] A5007 line. Furthermore, the FWHMs 
of our synthesized emission-line profiles are significantly 
broader than 140 km s" 1 . These two discrepancies dis- 
favor the interpretation of CXOJ033831.8-352604 as the 
tidal disruption of a WD by a nd IMBH. The ob serva- 
tions of optical emission lines by IZepf et all (|2008| ) , how- 
ever, are consistent with the model in terms of both lu- 
minosity and FWHM. The observed luminosity of the 
[O in] A5007 line falls within the range of luminosities 
predicted f or system s with IMBH masses between 100 - 
1000 M Q . |Figure 9| shows a synthesized [O in] AA4959, 
5007 emission-line profile that has been smoo thed to the 
same s pectral resolution as the observations of Zcpf et al. 
(20081). The line profile was synthesized for model C with 
the observer at io = 89° and 9 = —12°. We have chosen 
these values to match the shape of the observed emission- 
line profile. The first similarity between the model and 
the observations is that the lines are so broad that they 
are blended. The FWHM of the model line profile is 
1700 km s _1 . Also, both the observed and modeled pro- 
files are asymmetric, with the blue wing much broader 



than the red. This suggests that we are oriented with 
the flow such that t> ma x is nearly parallel to our line of 
sight. While the synthesized and observed profiles are 
qualitatively similar, spectra of high signal-to-noise ratio 
and higher spectral resolution are required for a more 
detailed comparison to determine whether or not the ob- 
served line profiles are consistent with the model. 

One major drawback to this interpretation of the data 
involves the duration of the X-ray flare and the ob- 
served [O in] A5007 luminosity. An X-ray luminosity 
of 8.5 x 10 39 erg s" 1 was measured for this source with 
ROSAT in 1992 ([Colbert fc PtaJdl2002l ). The change in 
luminosity between t he ROSAT observation a nd the ob- 
servation reported in iMaccarone et al.l ((2007) is consis- 
tent with the model for the X-ray flare, if the ROSAT 
observation detected the flare ~ 1 year after the WD was 
tidally disrupted. However, this means that the [O m] 
A5007 luminosity was still 1.4 x 10 37 erg s -1 more than 
15 years after the tidal disruption, and long after the 
pea k in [O m ] A5007 emission predicted by this model 
(see |Figure 6[ ). The model can only explain this [O in] 
A5007 luminosity if the IMBH in this globular cluster is 
100 Mq and maximally spinning. Even though it is pos- 
sible to tune the model to explain the observed X-ray 
and [O in] A5007 luminosities, because the high emis- 
sion line luminosity persists for many years longer than it 
would in all but the most favorable scenario, it is unlikely 
that this source is a WD that has been tidally disrupted 
by a IMBH, according to the model. Furthermore, the 
la rge amplitude X-ray var iability of the source reported 
in lMaccarone et aLl (l2010f) is also not consistent with the 
model and lSteele et al.l (|2010f) find that the profile of the 
emission lines in the [O m] AA4959, 5007 doublet are 
well fit by a two component model with a bipolar conical 
outflow and a lower velocity Gaussian component. 

If, in spite of the above discrepancies, this source is a 
WD that has been tidally disrupted by a IMBH, then 
the strongest line in a UV spectrum will be C iv A1549. 
Unlike the [O m] A5007 emission line, the C iv A1549 
flux is fairly constant across the illuminated face of the 
cloud so most of this line's luminosity is generated in the 
outer portions of the debris tail where the emitting area 
is largest. This results in a emission-line profile consist- 
ing of a redshifted core with a FWHM of several hundred 
km s -1 and a broad blue shoulder. Furthermore, if the 
X-rays decline dramatically in the future, the [O in] lu- 
minosity should also dim over the course of one year ac- 
cording to this model. The delay in the decline is due to 
the light travel time from the IMBH to the distant por- 
tions of the unbound debris that generate most of the 
[O in] luminosity. 

Given the low rate at which WDs are tidally disrupted 
by IMBHs, the short time during which the disruption is 
observable, and the simplifying assumptions made in the 
model presented here, it is not surprising that the emis- 
sion from CXOJ033831. 8-352604 and the globular cluster 
associated with NGC 4472 are not completely consistent 
with the predictions of the model. On average, detect- 
ing a single WD tidal disruption in a globular cluster 
requires surveying 10 6 galaxies, so it is unlikely that one 
has been observed. Additionally, improved models that 
consider the dynamics of the unbound debris in more 
detail and include emission from the material ejected in 
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TABLE 1 
Model Parameters 







Rp 


Rlso 


*Edd 


A</) h 


L ma x 


Model 


(Mq) 


R T 


Rs 


(days) 


(rad) 


C 


A 


10 2 


1 


3 


1100 


1.5 


0.032 


13 


10 2 


0.3 


3 


240 


2.7 


0.083 


C 


10 3 


1 


3 


440 


0.99 


0.052 


D 


10 3 


0.3 


3 


95 


1.9 


0.15 


E 


10 3 


1 


0.5 


440 


0.99 


0.052 


F 


10 4 


1 


3 


170 


0.68 


0.082 



a The length of the supcr-Eddington mass fallback phase 
b The azimuthal dispersion of the unbound debris in the 
orbital plane 

c The maximum velocity in the unbound debris cloud as 
a fraction of the speed of light 



the super-Eddington outflows could capture behavior in 
the optical emission line spectrum that is not predicted 
by this model. However, given the capabilities of cur- 
rent and planed optical and X-ray transient surveys, the 
number of candidate WD tidal disruption events should 
increase in the near future. If follow up optical and UV 
spectroscopy do not reveal strong, broad [ O hi] A50 07 
and C iv A 1549 lines that evolve as shown if |Figure~6} it 



is unlikely that the emission is produced by a WD that 
has been tidally disrupted by an IMBH. 
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TABLE 2 
Predicted Line Luminosities 



Luminosity of Luminosities of Selected Lines Relative to [O III] A5007 

[O III] A5007 

Model (ergs" 1 ) [O III] A4363 [O II] A7325 C III A977 C IVA1549 [C I] A8727 

?icr (cm" 3 ) 2 -► 7 x 10 5 3 x 10 7 6 x 10 6 ~ 10 7 



600 days 



A 


1.2xl0 37 


0.094 


0.17 


0.12 


3.4 


8.3xl0" 3 


B 


1.1X10 36 


0.37 


0.07 


0.11 


6.6 


9.4xl0~ 5 


C 


5.6xl0 36 


0.15 


0.16 


0.22 


8.1 


3.7xl0~ 3 


D 


5.9xl0 35 


0.044 


0.055 


0.17 


19 


4.0xl0~ 4 


E 


2.5xl0 37 


0.12 


0.33 


0.14 


3.9 


1.9xl0~ 2 


F 


2.7xl0 36 


0.17 


0.15 


0.17 


5.9 


5.6xl0~ 3 


2500 days 


A 


3.9xl0 36 


0.018 


0.035 


0.059 


2.2 


3.1xl0" 3 


B 


1.2X10 38 


0.010 


6.7xl0~ 3 


0.077 


6.9 


6.3xl0~ 5 


C 


2.1xl0 36 


0.016 


0.029 


0.054 


2.0 


2.9xl0~ 3 


D 


4.3 xlO 34 


0.013 


0.014 


0.16 


26 


3.8xl0~ 4 


E 


9.0xl0 36 


0.019 


0.071 


0.055 


1.9 


5.1xl0~ 3 


F 


8.6xl0 35 


0.015 


0.044 


0.051 


0.80 


3.6xl0" 3 



a The crit ical densities of the forbidden tr ansitions were taken from the following sources: [O III] 
lines from IDe Robertis fc OsterbrocH iT986ft [O II] A7325 from ITsyetanov fc Iankuloval (fT989l1. The 
value for [C I] A8727 was es timated from the data provided by INussbaumer fc Ruscal l|1979fl and 
IPeauignot fc Aldrovandil U976! ). 
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-0.10 -0.08 -0.06 -0.04 -0.02 0.00 0.02 

x [pc] 

Fig. 1. — The spatial distribution of the unbound debris from model C, projected into the orbital plane. The BH sits at the origin and the 
WD approached it from the — x direction on a parabolic orbit. The debris is shown one, three, and seven years after the tidal disruption. 
The shading shows the density of the debris, which goes at R~ 3 . The ^-curve begins at the BH and lies along the illuminated face of the 
debris cloud. The ionizing radiation from the disk encounters the cloud at a different angle of incidence for each point along i and travels 
along the corresponding h into the cloud. The length of h shown here is greatly exaggerated, the incident radiation only reaches depths 
< 10~ 3 R. 
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X [A] 

1 5 1 4 1 3 1 2 1 1 1 0° 1 5 1 4 1 3 1 2 1 1 1 0° 




hz/ [keV] 

Fig. 2. — The parameterized AGN (solid line), multicolor blackbody (dashed line), and the multicolor blackbody plus X-ray power law 
(dotted line) SEDs incident on an azimuthal segment of the debris cloud considered in model C. The segment has R = 1.9 X 10 16 cm and 
R = 1.6 X 10 17 cm at t = 300 days and t = 2500 days, respectively. Note that over time the position of the peak moves towards longer 
wavelengths as the temperature in the disk decreases and the total flux decreases as the accretion rate declines and the segment moves 
further away from the ionizing source. 
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10 15 10 16 10 17 
t [cm] 

Fig. 3. — The oxygen ionization parameter Uo (dashe d line), ra dial distance from the IMBH R (dotted line), and the number density n 
(solid line) along £, the illuminated face of the cloud (see |Figure 1) . These quantities are plotted at 600 (blue) and 1200 (red) days after the 
tidal disruption for model C. At these times, the super-Eddington outflows have subsided so Uo declines as t~ 0,2 . Note that R is nearly 
proportional to I. The horizontal, dot-dashed line shows the critical density of [O III] A5007. At later times, a larger portion of the cloud is 
below n cr and the region of maximum emissivity, the portion of the cloud with n = n cr , shifts to a part of the cloud moving with a lower 
velocity with respect to the IMBH. 
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Fig. 4. — The fractions of the first four ions of oxygen (left panel) and carbon (right panel) along h. These values are from a photoionization 
calculation for an azimuthal segment near the middle of the debris arc with R = 1.8 X 10 17 cm, n = 7.5 X 10 4 cm -3 , and Uo = 0.08, 1500 
days after tidal disruption in model C, using the multicolor blackbody plus X-ray power law SED. Highly ionized species are present in the 
surface layers of the cloud and singly ionized carbon and oxygen are the dominant ions in deeper layers. In most of our photoionization 
models the dominant ions of oxygen and carbon are O II and C II. 
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Fig. 5. — The light curves of five of the strongest emission lines. These light curves come from dynamical model C and were calculated 
using the multicolor blackbody plus X-ray power law SED. The dotted portion of each curve shows the line luminosity while the fallback 
rate is super-Eddington and our model is uncertain. The dashed line illustrates the t -5 / 3 dependance of the mass fallback rate. The 
normalization of the dashed line is arbitrary and is shown for comparison with the slopes of the emission line light curves. 
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100 1000 100 1000 

time [doys] 

Fig. 6. — The [O III] A5007 (left panel) and C IV A1549 (right panel) light curves for each model given in ITable T1 The photoionization 
calculations that produced these light curves used the multicolor blackbody plus X-ray power law SED; calculations that exclude the 
X-ray component produce nearly identical results. The dotted portion of each curve shows the line luminosity while the fallback rate is 
super-Eddington and our model is uncertain. The curves in the left panel differ from those in the right because [O III] A5007 is a forbidden 
emission line and the density in the debris cloud must drop to the critical density for this transition for this emission line to form. Curves 
B and D in the left panel are different from the others because they correspond to models with R p = 0.3Rt- Similarly, the curves in the 
right panel corresponding to models with R p = 0.3Rt differ from all the others by their steep decline at early times. The dashed line 
illustrates the 4 -5 / 3 dependance of the mass fallback rate. The normalization of the dashed line is arbitrary and is shown for comparison 
with the slopes of the emission line light curves. 



EM Emission From Tidally Disrupted WD 



17 



in 
-t—i 

"c 

D 

o 

■A— I 

i_ 
o 



CD 

> 
o 

CD 

cr 




0.01 0.02 0.03 0.04 



Fig. 7. — The time evolution of the [O III] A5007 emission-line profile as seen by an observer at iq = 45° and 9q = —5°. Time increases 
from bottom to top and the time for each profile is written on the plot. The dotted lines show the = level for each profile. The 
horizontal axis shows z = (A — 5007A)/5007A. Each profile has been normalized so that it has a maximum height of 1 and the profiles are 
offset by one unit for clarity. Note that in addition to becoming broader with time, the peak of the line profile shifts to lower velocity. Once 
most of the cloud is below n cr , the emissivity per unit volume decreases throughout the cloud and the peak shifts back towards higher 
velocities that correspond to regions in the tail with larger volume and a larger fraction of O III. 
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Fig. 8. — The time evolution of the [O III] A5007 emission-line profile as seen by an observer at iq = 45° and 6q = 30°. Time increases 
from bottom to top and the time for each profile is written on the plot. The dotted lines show the = level for each profile. The 
horizontal axis shows z = (A — 5007A)/5007A. Each profile has been normalized so that it has a maximum height of 1 and the profiles are 
offset by one unit for clarity. Note that in addition to becoming broader with time, the peak of the line profile shifts to lower velocity. Once 
most of the cloud is below n cr , the emissivity per unit volume decreases throughout the cloud and the peak shifts back towards higher 
velocities that correspond to regions in the tail with larger volume and a larger fraction of O III. 
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Fig. 9. — The [O III] AA4959, 5007 line profile for model C, 1750 days after tidal disruption, as seen by an o bserver at = 8 9° and 
00 = —12°. The profile has been smoothed to a spectral resolution of 400 for comparison with observations of IZepf et al, (2003). The 
profile is normalized to have a height of 1 and does not include contributions from the continuum. 



